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The center-bound excitonic diffusion on dendrimers subjected to several types of non- 
homogeneous funnehng potentials, is considered. We first study the mean-first passage time (MFPT) 
for diffusion in a linear potential with different types of correlated and uncorrelated random per- 
turbations. Increasing the funneling force, there is a transition from a phase in which the MFPT 
grows exponentially with the number of generations g, to one in which it does so linearly. Overall 
the disorder slows down the diffusion, but the effect is much more pronounced in the exponential 
compared to the linear phase. When the disorder gives rise to uncorrelated random forces there is, 
in addition, a transition as the temperature T is lowered. This is a transition from a high-T regime 
in which all paths contribute to the MFPT to a low-T regime in which only a few of them do. 
We further explore the funneling within a realistic non-linear potential for extended dendrimers in 
which the dependence of the lowest excitonic energy level on the segment length was derived using 
the Time-Dependent Hatree-Fock approximation. Under this potential the MFPT grows initially 
linearly with g but crosses-over, beyond a molecular-specific and T-dependent optimal size, to an 
exponential increase. Finally we consider geometrical disorder in the form of a small concentration 
of long connections as in the small world model. Beyond a critical concentration of connections 
the MFPT decreases significantly and it changes to a power-law or to a logarithmic scaling with g, 
depending on the strength of the funneling force. 



I. INTRODUCTION 



Extended dendrimers (Fig 1.) are nanoscale Cayley tree like supermolecules which exhibit an energy gradient from 
the periphery to the center of the molecule. Their unique hierarchical self-similar structure can be described by three 
ingredients: the basic building block (e.g. phenyl acetylene linear segment), the branching of the end points and the 
number of generations. The number of basic elements grows exponentially with the number of generations g and the 
number of elements at the periphery of the molecule is comparable to the number of bulk elements. 

Recent theoretical and experimental studies |^ ^ Q| have shown that the electronic excitations of these den- 
drimers are spatially localized within each segment. For these dendrimers, lengthening of the linear tree branches 
towards the core leads to a hierarchy of localization lengths. Hence the exciton energy decreases with generation 
from the periphery to the core. This energy funnel combined with the large number of absorbing elements at the 
periphery make these dendrimers potentially efficient single molecule antenna systems Excitons created upon 
optical excitation of the shortest linear segments at the periphery, diffuse through the intermediate regions, finally 
reaching the core where an energy trap is located. Excitation transfer proceeds via Coulomb interaction and may be 
described by the Frenkel exciton model. Theoretical studies focused on calculating the time it takes for an exciton 
generated at the surface to reach the core. This trapping time is a direct measure of the efficiency of the antenna. Its 
dependence on the number of generations g and the funneling driving force were calculated. 

These studies which were able to capture some important features of excitonic diffusion on dendrimers, assumed 
that all the tree branches (namely the linear polymeric chains) of a given generation are identical. This implies 
that the exciton energies are fixed. In reality interaction with the solvent and intramolecular vibrations induce slow 
(quenched) fluctuations in the energy e. (Nonlinear Q and single molecule spectroscopy typically show nanosecond 
to millisecond bath motions whereas the exciton trapping times are typically in the picosecond range pi | p|). 

On a more fundamental level, diffusion in disordered media is an active field of Statistical Physics |17t. Different 
forms of correlations among energy levels were considered and nontrivial behaviors, such as the difference between the 
typical and the average diffusion, anomalous scaling, and others were found in some cases. Independently, theoretical 
investigations of dynamics on Cayley trees (or Bethe lattices) lead to many interesting results. Exact solutions exist for 
some problems and bear important consequences. Directed polymer on Cayley trees is one of the interesting examples 
where a spin-glass like transition was predicted pO{ . As shown below the problem of diffusion on a disordered Cayley 
tree exhibits such a transition as well, and thus has a considerable theoretical significance. 

In a recent letter [ pT| , we reported preliminary studies of the effects of nonlinear dependence of the excitonic energy 
and realistic quenched disorder on excitonic diffusion on Cayley tree like dendrimers. We further carried out analytic 
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investigations of diffusion on a disordered Cayley tree with a linear potential. Our findings can be summarized as 
follows: 

(i) Due to the non-linear variation of the funneling potential with the generation, an optimal the generation number 
n*{g) was found beyond which the nonlinear potential drastically diminishes the light-harvesting efficiency. Moreover, 
increasing of the trapping time by disorder (by slowing down the exciton diffusion) will be more pronounced if the 
dendrimer crosses this optimal size. Hence, for dendrimers larger than the optimal size, even though the total photon 
absorbance will increase due to the increase in the number of peripheral light absorbing sites, the slower excitonic 
migration towards the active center will make the light-harvesting antenna less effective. 

(ii) We investigated the diffusion on a Cayley tree in the presence of a linear potential and considered various types 
of disorder, depending on the correlations among the energy fluctuations. For a specific type of disorder (See model 
(iv) below), we found a dynamic phase transition to a highly disordered phase where only a few paths dominate 
the exciton migration. This resembles the equilibrium replica symmetry breaking transition found in other random 
systems (e.g. directed polymers in random media). 

The purpose of this paper is to elaborate on, and extend the results reported in our previous study. We further 
report the numerical study of diffusion on a Cayley tree with a specific kind of geometric disorder obtained by adding a 
few random connections between various sites to the usual tree structure. This small world model was first introduced 
by Watts and Strogatz for Euclidean lattices. 

The organization of the paper is as follows: The MFPT for random walk is reviewed in Sec. II. In Sec. Ill we 
discuss diffusion on a Random Cayley tree with a linear potential. The effect of nonlinear potential on the exciton 
migration is reported in Sec. IV. The small world model is studied in Sec. V. Finally, we conclude in Sec. VI with a 
summary of our results. 



II. RANDOM WALKS AND THE MEAN FIRST PASSAGE TIME (MFPT) 

The Mean First Passage Time (MFPT) is a faithful measure of the trapping efficiency defined as the average time 
it takes for an exciton to diffuse from the periphery to the core, where it gets trapped [|l3|. Another quantity of 
interest related to the MFPT is the Mean Residence Time (MRT), which is the average time spent on a site of the 
tree (branch of the dendrimer)[13,14,23,24]. Our theoretical effort focused on calculating the MFPT and the MRT 
for a continuous time random walk with exponential distribution of waiting times which gives a master equation for 
the probability Pn{t) of the exciton to be on site n at time t. If the energy depends only on n, the problem becomes 
effectively one-dimensional for Pn{t). The random walker on the n-th site of the one-dimensional chain jumps to 
nearest sites with rates Ri (toward the reflecting point) and Ti (toward the trap). The random walk is described by 
the master equation |l^] 

where P{t) denotes the probability vector and the survival probability S{t) = J2i Pn{t)- A is the tridiagonal transition 
matrix. Written explicitly in terms of the rates R and T: 

Po{t)=T,P,{t), 

Pi{t)^T2P2{t)-{T^ + Ri)Pi{t), 

Pn{t) = Tn+lPn+l{t) + Rn-lPn-l{t) " {T^ + Rn)Pn{t) (K n < ff), (2) 

Pg{t)=Ra-lP9-l{T)-TgPg{t)- 

The trap (Pq) is not included in Eq. (1), while reflecting Boundary Condition are imposed &i n = g. 
The formal solution of Eq. (1) is 

V{t) = exp(-A<)P(0), (3) 

where P(0) is the initial condition. P(0) = [0, 0, 0....1, 0, 0, 0...]"'" for an excitation which starts on the n-th generation. 
The MFPT is obtained by integrating the survival probabilities P„(t) over time 
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where tn is the MRT for generation n. Using the general solution Eq. (3) [ p^ 

/ P{t)dt^ / exp{~At)P{0)dt = A-^P{0). (5) 
Jo Jo 

Hence the MRT at site n is given by 

/•oo 

tn^ p„(t) = 5]4;/p,(0). (6) 

Jo ^ 

For an excitation which starts at the highest (periphery) generation g, Pj{0) — Sjg and the MRT is given by 

ti = ^ n = 1 (7) 

^ n—2 ^ n— 1 

^« = ^ + E(r^ n l<n<g (8) 

" 1=0 * j=i+l 

where ^„ = i?„/T„+i is the detailed balance ratio. The MFPT is given by 

9-1 

r{g) = E (9) 

9 , 9-2 m-2 m-1 

= E^ + EE(^) n (10) 

m=l m=2 2—0 

An explicit form of this equation is given in the appendix. Hereafter we assume r„ = 1 and treat ^„ as independent 
variables with ^„ = cexp{— /3(e(n + 1) — e{n))} where c + 1 is the coordination number of the tree (c = 2 in our case), 
e{n) denotes the energy of the nth segment and (3 = (KT)^^ is the Boltzmann factor. 

III. RANDOM WALK AND MFPT ON A DISORDERED CAYLEY TREE WITH A LINEAR FUNNEL 

The linear tunneling potential assumes the form 

e{n) = nF (11) 

F being the potential difference. For this model ^„ = Co = cexp(— /3F). Eq. (10) can be easily summed (see the 
appendix for details) to yield for the MFPT 

For large g (^ 1) Eq. (12) shows three distinct behaviors 

(a) C > 1: T{g) ~ ex.p{ghi£^) (exponential regime). 

(b) < 1- '''(5) ^ 9 (linear regime). 

(c) C = 1: T{g) ~ g^ (diffusive behavior). 

Real dendrimers have energy fluctuations and quenched (slow) energy disorder plays an important role in the 
energy transfer dynamics. We have studied the effects of disorder on diffusion in the presence of a linear potential by 
considering four models of disorder denoted by (i) - (iv). These will be described below. 

A. Intergenerational Quenched Disorder 

Assuming that energy fluctuations in the segment lengths of the same generation are identical, the diffusion is 
mapped into an effective one-dimensional problem. We denote this type of disorder intergenerational. 
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Intergenerational energy fluctuations can be introduced in two ways. In the absence of correlations, we obtain the 
standard diagonal disorder (random energy) model (i) where a random part e„ is added to the linear energy e(n) . In 
the second form of disorder (random force model(ii)), the energy differences Ae„ (= e{n + 1) — e{n))) are treated as 
random variables with the positive average values e(n) to ensure tunneling. Both models are one dimensional and the 
MFPT is given by eq.(lO). 

Model (i): Random intergenerational energy 

In this model, all segments in a given generation have the same random energy. Generational energies are made 
random by adding a fluctuating part e„ to the linear potential 

e{n) = e(n) + e„. 

e„'s are independent and identically distributed with the average (e„) = and e{n) is the linear potential defined in 
Eq. (11). A similar random energy model has been considered previously ||l^ for a discrete random walk. 

Hereafter () will denote the average over disorder. We define r]n — exp(— /3e„) and assume {r]^^) and (77^^) to be 
finite. The subscript n in (77) is omitted because the random variable e„'s are independent and identically distributed. 
Hence 

= cexp[-p{e{n + 1) - e{n))] cexp[-(3{e{n) - e{n - 1))].... (13) 

...cexp[-p{e{3) - e(2))] cexp[-P{e{2) - e(l))] 
= ^^{exp[-pien+i - £„)]}. (14) 

Averaging over realizations of disorder gives 

< e.e.-i.-6a >=< V X > (15) 

Substituting this in Eq. (10) gives for the MFPT 

(rig)) V >< V'^ > Ta{g), (16) 

where To{g) is the MFPT (see Eq. (12)) for a linear potential with no disorder. 

Similar to the ordered case, the disorder averaged MFPT has three distinct regimes: linear (^0 < 1), quadratic 
(^0 = 1) and exponential (^0 > 1) depending on the value of ^q. Even though the g dependence and the critical 
point (^0 = 1) of the MFPT do not change by disorder, the magnitude of the disorder averaged MFPT exceeds the 
corresponding value for the ordered system. Disorder slows down the first passage diffusion. If the fluctuating part 
of the energy e„ has a Gaussian probability distribution Pg{() = (l/\/27rA)ea;p{— e^/2A^}, then 

(Tg) - < 77 >< 77-1 > ro(<7) 

= {exp{-f3en)){exp{(3en))To{g) 

= exp{(3^X^)To{g). (17) 
The average MFPT is thus increased by the disorder-dependent factor of exp{f3'^ X^) . 

So far we have considered only the average MFPT, but due to disorder there will be fluctuations in the MFPT 
corresponding to various realizations of disorder. The distribution of MFPT due to disorder may be directly observed in 
single molecule spectroscopy |l^. We have calculated the rms fluctuations in the MFPT < (Ar)^ >=< (t— < t >)^ > 
in different regimes of < r(g) >. 

In the linear regime (^0 < 1), only the terms linear in ^ in Eq. (10) arc sufficient to yield the correct scaling (in 
terms of g) for both < T(g) > and < Ar^ >. Considering only the linear terms in Eq.(lO), we get 

rig) ~ {1} + {Co exph/5(e2 - ei)] + 1} + Uo exp[-/3(e3 - £2)] + 1} + .... 
.... + {eoexph/3(eg-eg_i)] + l} 
= ^o{exp(-/?e2) exp(/3ei) + exp(-/3e3) exp(/3e2) + .... 

+ exp(-/3eg) exp(/3eg_i)} + .g. (18) 

Averaging over realizations of disorder gives 

{r{g)) = {g-i)io{v){v-')+g (19) 
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To calculate the fluctuation of the MFPT around its average we need ((r((7))^). Taking the disorder average of the 
square of the Equation 

{(rigyf) = [g - 1) e, ivfiv-') + 2(3 - 1) eo iv) iV') 

+Mg-'^)^o{v){v-')+g'- (20) 

From the average MFPT eq. (19) 

{{rigW - (.9 - Ifeo ivriv-y 2.g(.g- 1)^0 (ri)(rr')+g'- (21) 
Combining eqs. (20) and (21) we obtain for the fluctuation in the MFPT 

((Ar(g))2) = ((r(<?))^) - {(rig)))' 

^geo {{v'){v-')+2{v){v-')-3{vr{v-r}- (22) 
Hence the relative fluctuation in the linear regime (^o < 1) will be given by 

At ((Ar(5))2)V2 



(rig)) 

g{i + Uv){v-')} 



(23) 



At the transition point = li we have to keep all terms in Eq. (10) in the calculation of Ar/r. Some tedious 
calculations finally yield 6t/t ~ l/y^, the same scaling as found above in the linear regime. 

In the exponential regime (^o > l); the MFPT is dominated by the single term in expression (10). In this case 

T{g) - C9-iCs-2--.66 

The average MFPT is then < T{g) >~ (,q < V >< -q^^ > for large g, as expected. In this approximation 
< (T(g))2 >~ < T]^ X 77^2 >, and the fluctuation in the MFPT is 



At r< ?7 >< T]^ > 



T ' (< ?7 >< 1]-'^)^ 



1]2 (24) 



which is a constant independent of g. 
Model (ii): Random intergenerational force 



In this model energy-differences between consecutive generations are randomly distributed, but again all segments 
in the same generation are identical. Fluctuations of the energy levels are correlated and the energy differences Ae„ 
are assumed to be identically distributed according to some probability distribution P(Aen). This diffusion model 



was studied extensively in various contexts and was reviewed in |20|. The effect of disorder is more pronounced in 



this case compared with the random energy model. The average MFPT is given by 

< r(g) > again has three regimes, but the transition occurs at < ^ >= 1. Note that in model (i), this condition 
was ^0 = Ij a-iid disorder did not change the transition point. In the present model, however, the transition point is 
shifted due to the stronger effect of disorder. In some regimes of the MFPT, the effect of disorder is so strong that 
the typical diffusion time differs from the MFPT. This indicates a broad distribution of the MFPT where the average 
is affected by the rare configurations with long first passage time. The more representative "typical" diffusion time is 



calculated as Ttyp = exp{logT{g)) . Following |20|, we define different scaling regimes of the MFPT and Ttyp: 

(a) For < ^ > < 1, < r(g) > is linear in g. <^ > determines the fiuctuation around the average. < T'^{g) >~ g^ as 

long as < >< 1 and hence the relative fiuctuation St/t goes to a constant (similar to the random energy case). For 

< >> 1, however, < T^{g) > grows exponentially with g as the typical behavior starts to differ from the MFPT. 
If < ^ > > 1, < r(g) > is always exponential in g. However, Ttyp varies depending on the value of < In^ > (whether 

(,=,or)l). 
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(b) < ^ > > 1, < /nf > < 0: The typical first passage time {Ttyp{g) = exp < logT{g) >) does not follow the 
average exponential behavior, but rather grows like a power law g°' with a =< S(ln^)'^ > /2 < Zn^ >. For a Gaussian 
probability distribution PciAe) the above formula yields a = fi^X"^ /2{ln^o). 

(c) < ^ > > 1, < ln£, >= At this transition point (the "Sinai point"), < T{g) > is exponential in g, but the 
more representative Ttyp{g) is only exponential in y/g [ p^ . 

(d) < <^ > > 1, < In^ > > 0: In this regime, both < T{g) > and Ttyp{g) diverge exponentially with g. 



B. Intersegment Quenched Disorder 

When energy fluctuations exist within the same generation (intersegment disorder), the system can no longer be 
mapped into a one dimensional model and we have to use eqs. (4) and (6) to compute the MFPT for the actual 
tree structure. The energy of site qn is denoted eg(n). The random variables ^„'s become matrices f^''? where {p,q) 
represents gth point of nth and thepth point of the (n + l)th generation, (note that c is not included in this definition). 
The detailed balance ratio becomes ^^'^ = cxp{— /3Aepg(r7,)} = ^ exp{—f3{ep^_^_^ — e^^)} , where eg„ is the fluctuating 
part of the energy and ^ = exp{— /3(e(ri + 1) — {e{n))} (so — c^)- From now on we assume the branching ratio 
c = 2. The MFPT for a particle released at the peripheral site Ig of the g-generational tree is 

(r'^(5)>={l} 



+ [e,' {e/ {■■■■{)■■■) +e/{ o....)) (^^'^^(....o....) +r ( o-.)) + •• 

+ (Cg-i+eg-i) + l} (26) 



We define the initial site average MFPT (ISA-MFPT) as: 



3 ig=l 

where Ng = [1 + cjc^^^ is the number of peripheral sites (leaves of the tree). 
T{g) may be expressed again (see Eq. (9)) as a sum over ISA-MRT: 

71=1 

Each of the tn is made of sum of products over ^'s up to generation n. Formally it may be expressed as 

^=i + E{^E^a' (29) 

fc=l ifc=l 

where T^" stands for the contribution to i„ from a subtree rooted at site ik on generation k and spreading out to the 
^th generation. This contribution is made of products of {n — k) ^'s along the c""'^ paths going from the site ik to 

the c""*^' different sites on generation n. The fc'th generation has c^~^ rooted subtrees with = 1, ,c'^~^. The 

contributions from subtrees rooted at the same generation are independent from each other. Symbolically, 

^".=e( n 

= ^exp -/3^Aep,(m)|, (30) 
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where Tm are all paths on this subtree from ik to generation n. 
Model (iii); Random intersegment energy 

In a fully disordered tree all segment energies (ej„)are random and uncorrelated . The techniques applied in model 

(i) can be simply extended to this case with almost identical results (Fig. 2 shows the average MFPT for both the 
models). 

In the linear regime, the MFPT expression can be approximated by neglecting terms and higher. For c = 2 we 
have 

r'^{g)^{l} 

+ exp[-/3(ei, - eij] + ^ exp[-/3(ei, - ei J] + 1} 
+ exp[-/3(ei3 - eij] + ^ exp[-/?(e23 - eij] + 1} 



+ exp[-/3(ei, - ei^_J] +e cxp[-/3(e2, - ei,_J] + l} 

= ^{exp[-/3(ei2 - eij] + exp[-f3{ei^ - eij] + exp[-/3(ei3 - eij] + exp[-/3(e23 - eij] 

+ + exp[-/J(ei^ -ei^_J]+exp[-/3(e2, -ei^_J]}+5. (31) 

Upon averaging over realizations of disorder we obtain 

{T'^{9))=2{g-l)^{r^){n-')+g 

= (5 -1)^0 +5 for g large, (32) 

where — 2C with c = 2. 

The average MFPT is identical for models (i) and (iii). The fluctuation of the MFPT around its average are 
connected to {{T{g)) ). Taking the disorder average of the square of the Eq. (31) 

{{r{g)f) = 2(5 - 1) f iv'Hv-') + 2(5 - 1) e ivfiv-') + 2{2(ff - 1)} e {v){v-') 

+2 { - - 1) - - 2)} e ivfiv-r 

+2{2gig-l)}^{r^){r^-')+g\ (33) 
The average MFPT eq. (32) yields 

{iT{g))f = 4(5 - 1)^ e - 4ff(5 - 1) e {ri) {r') + 9^ 

= (5 - 1)' Co - '^gig - 1) ^0 {ri){v-') + 9' (34) 

Eqs. (33) and (34) yield for the fluctuation in the MFPT 
((Ar(5))') = {{T{g)f) - {{T{g))f 

= g^o [l{v''){v-') + livfiv-') + {v){v-') - 2{vf{v-r] (35) 

The relative fluctuation Ar/r ~ l/y^- Scales with g the same way as in the case of intergenerational disorder, but 
the amplitude is smaller in this model. This can be; seen by comparing Eq. (22) and Eq. (35) and since {r])^ < {rf). 
This is also verified by our numerical calculation (Fig. (3)) where we kept all the terms in the A^IFPT expression. In 
the exponential regime, the relative fluctuations At/t will still saturate to a (^-independent value, but the latter will 
be also smaller than that reached with Model (i). 

The differences between models (i) and (iii) may also stem from the additional fluctuations in the MFPT for the 
latter model which arise from distinct initial sites at the periphery. For initial excitations starting at two different 
peripheral sites (on the same dendrimer) separated by an ultrametric distance p the fluctuations in the MFPT is given 
by 

< (Ar(5))2 >~p. 

Averaging over all possible peripheral sites, the fluctuation scales with g the same way as the dendrimer-to-dendrimer 
fluctuations discussed above. The amplitude of the fluctuation is different from that of the dendrimer-to-dendrimer 
fluctuations and depends heavily on the branching ratio c of the tree. 
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Model (iv): Random intersegment force 

This model assumes that all energy-differences (Aepq(n) = ep(n + 1) — e.q{n)) between neighboring segments are 
random. We observe a new transition apart from the usual linear to exponential regime transition of the MFPT. 
This is a one-step replica symmetry breaking transition from a weakly disordered (high-temperature) phase to highly 
disordered (low-temperature) phase. In the weakly disordered phase all paths contribute to T(g), whereas in the 
highly-disordered phase T{g) is dominated by a few paths. Similar dynamic transition was predicted ]2^ ] for kinetic 
pathways in protein folding. 

This transition is occurring only in this model. It cannot occur in models (i) and (ii) since in these models there 
is effectively a single pathway from the periphery to the center (i.e. all pathways are equivalent). It cannot occur in 
model (iii) since the random-energy effect is very limited. For any path, independent of its length, only the random 
energies at the initial and final sites are contributing. This is insufficient to cause an energetic disparity which will 
overcome the entropic advantage of having a maximum number of path contributing. 

In the low-temperature phase the distribution of — exp(— /3Aep ^(ri)) becomes very broad. The important 
contributions will come from those paths for which X](pq)„ ^^p,q{^) is large and negative. We also assume that the 
constant force e = is weak and the system is in the exponential regime. The linear regime will be discussed later. 
We thus focus on the large ^ regime (it is sufficient for them to be typically larger than 1/c). In this regime, the 
largest contributions will come from the longest paths, and if a qualitative change in the behavior will occur it will be 
noticeable first in the dominant contribution from the largest of all subtrees. We will thus look for a new behavior in 
the ISA-MKT by examining at the maximal tree rooted at fc = 1. The same argument may be repeated while looking 
at T{g): its dominant contribution in this regime will come from tg namely the ISA-MKT on the most extremal 
generation with n ~ g. 

Searching for a possible abrupt change in the properties of ISA-MFPT, we will therefore examine the largest tree 
contributing to tg ~ Tf since it has the largest energy disparity. As long as its behavior is normal, so will be all 
other terms. Once its behavior changes, it will affect that of tg and then that of T(g), to which it makes the largest 
contribution. 

The contribution of Tf is akin to a partition function of a so-called "random-directed polymer" on a Cayley tree 
which is known to have a glass transition (in the thermodynamic limit). Physically the transition is between a high 
temperature phase T > Tg at which all paths contribute, while for T < Tg only a finite number of them do. The 
transition is of one-step replica symmetry breaking (IKSB) type. The n — > replica-trick was used to study its 
transition point. A complete analysis using this trick is available JTo| . Here we present a simpler, more intuitive, 
outline. In the high temperature regime (3 < (3g \t may be shown that all moments of the MFPT (or the "partition 
function" , in the polymer picture) obey {T^) = {T)^ . Hence all thermodynamic properties of the quenched system 
are given by those in which the disorder is annealed. In this regime < Tf > is easily computed and found to be 

< Tf >= (c < exp(-/3Ae) >)». (36) 

For large g we have 

< Tf >= [exp{-/3/}]^ (37) 

where f{/3) = — In < Tj > is the free energy per generation. 

When the random energy difference Ae follows a Gaussian distribution Pc{Ae) = {l/^/27r\)exp{— ^^'^^l°^ }, the 
free energy is given by 

The transition point can be determined by the vanishing point of the entropy S{— — Inc — A^/2T^ , where we 
assumed Kb = 1). This yields Tg = X/\/2 logc, which is the transition temperature from the weakly disordered to 
the highly disordered phase. 

For low temperatures f3 > (3g the entropy remains zero. This is achieved ||l0| by expressing the free energy as 
f = eo — ^ — "'^'^ , and m g [0, 1] is chosen such that / is maximal. For P < Pg, m = 1 and the previous result 
holds. For P > f3g, §^ — (which implies ^ = 0) yields 



For a general distribution the requirement is to maximize 



/(to) = (Inc + ln(exp(— TO/3Ae))) . 

TOP 

The physical meaning of to is explained by the existence of a finite (in the thermodynamic 5 — s- oo limit) number 
of paths that contribute to Tf , and therefore to tg and T{g) as well. To see that, we need to define the overlap q 
between the contribution of two different products (paths) in the sum (and its distribution P{q)) which plays here 
the role of an order parameter, ft is defined as the weighted (over all paths) fraction of the way the two paths will go 
together on the same segments, averaged over the disorder. At high temperatures all paths are essentially equivalent 
and the probability for any overlap vanishes. At T = there is a single dominant path with the minimal energy (the 
"ground state" ) that makes the only contribution and thus both paths will overlap all the way giving q = 1 with 
probability 1. It was shown | [To[ | that for this system the only two possible overlaps are q = and q — 1. Their 
respective probabilities, however changes in the glassy phase and is given by: 

P(q) = 6{q) for P < Pg and 

P{q) = mS{q) + (1 - m)S{q - 1) for (3 > (3g (39) 

where m = S- was introduced before. 

So far the discussion was limited to the exponential regime. Increasing will lead to the transition into a linear 
regime at eg — For j3 < Pg the transition is given by the expected "annealed" condition c < ^ >= 1 which 

yields = ^ + \P^^ ■ For (3 > Pg, in the glassy phase free-energy is given by 

J./ m liic to/3A^ 

/to,/3=— + . 40 

Since to/3 = Pg the critical force remains constant within the glassy phase for Tg < T < and is equal to its value 
at the transition point (P = Pg), which is given by 



eoiP > /^.) - 4^ + ^ = A^/2h^^. (41) 
V 2 In c ^ 

The expected {eo,T) phase diagram is shown in Fig. (4) 

We should also note that crossing Pg within the linear regime yields only a glassy crossover and not a sharp 
transition. This is due to the fact that the ISA-MRT converges to a finite value even as g ^ 00. The longer paths 
thus make an exponentially small (in their length) contribution, the system is dominated by paths of finite length 
and the "glass transition" takes place in an effectively finite system. 

IV. MFPT FOR THE NON-LINEAR TDHF POTENTIAL 

In the previous section the funneling energy was assumed to vary linearly with the generation number n (= 0, 1, 2...g). 
However electronic structure calculations by Tretiak, Chernyak, and Mukamel show a strong nonlinear dependence 
of the exciton energy (on n) . 

The funneling effect of extended dendrimers originates from the variation of the segments length. We denote by 
n(= 1,2, ... ,g) the n*'* generation starting from the center. The peripheral n = g generation is made from m = 1 
monomer. The number of monomers increases by one going from one generation to the next towards to the center. 
The number of acetylene monomers in the n*'' generation is thus I — g ~ {n — 1). The n dependence of the excitation 
energy was computed using time-dependent Hartree Fock (TDHF) technique and fitted to the form |Q 

e(n) = A{1 + ^-^)°-' (42) 

(n- 1) 

with A = 2.80 ± 0.02{e.v.) and L = 0.669 ± 0.034. 
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A. Ordered Dendrimer 

Using eq. (42), the detailed balance ratio becomes 

«„ = c»p-A{(i,_i_)0.a.(, + _£_,o.a} 

The points of the same generation are identical and we can calculate the MRT and MFPT for this potential using 
the expressions (8) and (10) for one-dimensional random walk starting at the periphery. The mean time spent by an 
exciton in the nth generation is given by 

n 

K = exp[- A(l + J^j^f-'] E c("-^) exp[A(l + jz-j^t'] (44) 
The total time to reach the trap will be sum of all the t„ 's, and the MFPT is 

n=l n=l ^ ^ ' r=l ^ ^ ' 

Fig. 5 shows that the MFPT for this potential is quite different compared to that for the constant potential 
difference. We find that it depends linearly on g for the first few generations, but gradually changes to exponential 

with increasing g. 

The behavior of the MFPT can be better understood in terms of MRT, which is not a monotonic function of n (see 
Fig. 6) for the TDHF potential. At the generations near the periphery (g — n <^ g) the energy difference between 
generations is large and the funncling force is so strong that it overcomes the outward cntropic force of c > 2. Near 
the center, on the other hand, the larger is g, the weaker is the tunneling force and the entropic term dominates. The 
competing effects nearly cancel at some intermediate generation n = n* (g) where the exciton spends maximum time. 
For large g, an estimate for n* may obtained from the recursion relation satisfied (from Eq. (8)) by the MKT's 

tn = ^n-ltn-1 + 1 (46) 

For fjv to be maximal, t„+i/f„ =Cn+l/*n < l,hence^„ < 1. Similarly requiring t„/t„_i > 1 yields ^„_i-|-1/t„_i > 1, 
so for large r„_i (or equivalently large g) we may take > 1. r„ is therefore maximized as ^„ ^ 1 from below. 
Using expression (43) we obtain n ~ {g — k + 1) where k{= k{A, L, c)) is a constant (independent of n and g) and 
Ti*{g) is the closest integer. This formula works well for our system for g > 12 and holds approximately (off by only 
one generation) for 7 < g < 11. Clearly, the larger the g, the better is the estimate of n*{g). The reduced free energy 
u{n) = (3e{n) — nlogc is plotted in Fig. 7. The MRT is maximal at the site at which the free energy is minimum. 
For large g, both n*{g) and the energy difference between the center and n* increases with g. As a result, the time to 
reach the center from n* grows exponentially with g (the time to arrive at n* from n = g is always much shorter). If 
g is small enough, however, the value of n* and u{n*) increase only weakly with g and the MFPT dependence on g 
may be approximated by a power series (dominated initially by the linear term). 

For fixed g, the MFPT may be changed by variations of the coordination number (c) and the temperature (T). 
Decreasing the temperature makes the funneling more effective. As a result, the crossover from the linear to the 
exponential regime occurs at a higher n. Increasing c has a similar effect: both the outbound "entropic force" which 
competes with the energy funneling and the MFPT increase (Fig. 5). 



B. Random Energy Fluctuations 

We introduced disorder by adding a fluctuating part to the energy from an uniform probability distribution, similar 
to the random energy models («) and {in) for the linear potential. The maximum value of the fluctuation is taken to 
be ± 4% of the energy value. Intergenerational and intersegment types of disorder were considered separately. For 
the former case we used the MFPT expression (10) obtained for one-dimensional random walk in TDHF potential. 
For the latter, the MFPT formula (26) for a Cayley-tree can be used to calculate the time spent by an exciton in the 
n'th generation starting from a point (the points are labeled as 1, 2, 3...2^~^) at the periphery is given by 

n n—m 

i„ = E {exp[/?e(m, 1)] E exp[-/3e(n, r)]}, (47) 

m=l r=l 
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where e(n, r) is energy of the rth branch of the nth generation of the tree, r can take values 1, 2, 3 2'^~^. The MFPT is 

obtained by summing all the i„'s. The effect of disorder was studied by numerically averaging the MFPT over different 
realizations of the disorder. Fig. 8 shows disorder averaged (~ 10"* realizations) MFPT for both intergenerational 
and intersegment disorders and also the case with no fluctuation. Disorder seems to work against the tunneling and 
increases the MFPT. The effect of disorder is more pronounced in the exponential regime of the MFPT (Fig. 8). As 
for previously discussed disordered linear potential the disorder-averaged MFPT is similar for intergenerational and 
intersegment disorder, consistent with our analytical calculation for models (z) and (iii). However the fluctuations 
around the average MFPT ST{g), are different (Fig. 8) especially for large g and the relative fluctuation saturates to 
a larger value for intergenerational disorder. 

V. GEOMETRIC DISORDER: DIFFUSION ON A TREE WITH RANDOM CONNECTIONS 

Different forms of generic networks have drawn much recent interest. The small world model was introduced by 
Watts and Strogatz Q and received a considerable attention in the past two years |Q . The model assumes a regular 
lattice and fixed number of nearest neighbor connections with a few randomly chosen connections between vertices. 
The essential features of this model are: (i) High local connectivity which resembles a regular graph, (ii) The average 
distance between any two points scales as ~ logi, where L is the linear size of the graph. This is an important 
property of a random graph. 

Inspired by this model we have studied the effect of a few random connections on the exciton diffusion. We first 
considered the MFPT for a one-dimensional chain with some random connections made between any two vertices 
with a fixed small probability. For a given chain length L, there exists a critical probability for random connections 
{pc ~ L^/'' above which the small world effects are observed. We have performed a numerical simulation of the 

MFPT for a one-dimensional small world chain using a generalized form of the master equation 1., where we allow a 
few hopping between vertices chosen randomly. We then averaged over different realizations of random connections. 
When p > pc, we found that the average MFPT scales logarithmically (instead of linearly)with g (Fig. 9) for < 1 
. When ^0 > 1, the disorder averaged MFPT grows almost linearly instead of exponentially with g (Fig. 9). This 
may be rationalized in terms of the average distance between any two vertices which is supposed to behave as log 17 
, where g is the length of the system. Hence we can approximate a small world linear chain of length g hy a, regular 
linear chain of length logg. Then the average MFPT in a small world 1-dimensional chain will be linear in \ogg {i.e. 
T{g) ~ logy) for i^o < 1- For ^0 > Ij the MFPT is exponential in \ogg and hence T(g) ^ exp(alogg) = g", where a 
is a constant. In our simulations a ^ 1, so the MFPT grows linearly. As the probability p of random connections is 
decreased, the onset of the small world effect occurs at larger system sizes. 

On a tree we add a few random connections between vertices of different generations. The resulting diffusion on 
such a "small world tree" yields results similar to the one dimensional case. Disorder averaged MFPT again show 
logarithmic (^0 < 1) or power-law (^0 > 1) scaling with the number of generation g. 

VI. CONCLUSIONS 

We have investigated analytically the diffusion on a disordered Cayley tree in the presence of a linear potential. 
Fluctuations around the MFPT showed different forms of scaling (in terms of g) depending on the types of disorder 
and the scaling regime of MFPT. For a specific form of disorder (random intersegment force) we found a new dynamic 
transition in the MFPT. In the low temperature (highly disordered) phase, the MFPT is dominated by a few paths. 
This transition resembles the one-step replica symmetry breaking glass transition found in other disordered systems 
like random directed polymers. 

We have also considered exciton diffusion on dendrimers with a nonlinear funneling energy. For dendrimers with 
a nonlinear potential larger than a specific size, even though the number of light-absorbing sites at the periphery 
increases, the MFPT starts to grow exponentially (with g), resulting in slow exciton trapping. Quenched disorder slows 
down the exciton diffusion towards the center. This effect is more pronounced when the MFPT of the corresponding 
ordered system is in the exponential regime. Hence to achieve an efficient funneling of excitons, the number of 
generations must be restricted to some optimal value. We have determined this optimal size for a particular class of 
phenylacetylene dendrimers. For large dendrimers, the free energy attains its minimum at n*{g) where an exciton 
spends most of its time during its journey towards the center. Assuming a steady supply of long-lived photoexcitations, 
the excitons will start gathering at n*{g), rendering the single-exciton picture invalid. In that case, we have to consider 
exciton-exciton interaction and annihilation processes [25] 
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Finally, we have considered diffusion on a lattice with a few random connections ("small world") between any 
two vertices. The effect of random connections is to reduce the effective system size logarithmically. Numerical 
calculations show that the disorder averaged MFPT is also decreased by a logarithmic factor. Hence there will be 
an exponential gain in the cxciton trapping. This might have interesting consequences for information propagation 
through hierarchical social communication networks. 
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Appendix 

The MFPT Eq. (10) can be expressed as 
r{9) = {1} 

+ {a + 1} 

+ {66+6 + 1} 
+ 

+ + + 66-1 +6 + 1} 



{e3_i^3_2....66 + 6-iCs-2--6 + + 6-1 + 1} (48) 



where Tj = 1 is assumed. 

For the linear potential the sum in Eq. (10) assumes the form 

r(<?) = {l} 

+ {^ + 1} 

+ U'+C+1} 



+ {r +r"' + ••••+ c'+e+i} 

+ 

+ {^9-1+^9-2^ +^2+e + l} (49) 

(50) 

This can be recast as 

r{g) =g + {g-m + {g- 2)e + 2^-^ + (5i) 

This series can be summed to yield eq. (12) 
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FIG. 1. Extended Phenylacetylene Dendrimers with increasing total number of generations. 



FIG. 2. The MFPT for the linear potential with both types of random energy fluctuations (indistinguishable from each 
other) compared to that of the pure system (in the exponential regime). Inset: the same in the linear regime 



FIG. 3. The relative fluctuation in the MFPT for both types of random energy disorder (intergenerational and intersegment) 
with a linear potential (in the linear regime). 



FIG. 4. The {eo,T) phase diagram for the replica symmetry breaking transition 



FIG. 5. The MFPT vs g for the ordered TDHF nonlinear potential for different branching ratios of the tree 



FIG. 6. The MRT vs n for the ordered TDHF potential, for different values of g. 



FIG. 7. The free energy u{n) for the ordered TDHF nonlinear potential are plotted for different values of g. 



FIG. 8. The MFPT for the TDHF potential with both types of random energy (indiscernible from each other) compared to 
that of the pure system. {Inset: Their relative rms variations vs g). 



FIG. 9. The MFPT vs g for the small world model (exponential regime). Inset: same for the hnear regime. 
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